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The stability of a Nambu-Goto membrane at the equatorial plane of the Reissner-Nordstr0m-de 
Sitter spacetime is studied. The covariant perturbation formalism is applied to study the behavior 
of the perturbation of the membrane. The perturbation equation is solved numerically. It is shown 
that a membrane intersecting a charged black hole, including extremely charged one, is unstable 
and that the positive cosmological constant strengthens the instability. 
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qq . It is believed that in the early universe a series of vacuum phase transitions led to several types of topological 
£SJ ' defects p5| . Topological defects are relics of the early universe and are expected to convey some information on 
physics of very high energy scales beyond our reach with ground-based accelerators. On the other hand, topological 
defects are candidates for the seed of the observed large scale structure of the universe such as sheet-like or filamentary 
structures or voids. Thus topological defects are attractive examples which connect the high energy physics and 
cosmology and we might be able to have some information on high energy physics through cosmological observations 
such as gravitational wave detection in future. Topological defects, if they existed, would interact with other strong 
qq i gravitational sources such as black holes, and then they might have experienced large deformation and emitted some 
information on themselves as gravitational waves. If we succeed in detecting such gravitational waves and identifying 
them, we will be able to confirm the existence of topological defects and hence the occurrence of a vacuum phase 
■ transition. 

\ The topological defects with finite extent such as the cosmic string or the domain wall are known to produce 
O" 1 an unusual gravitational field fl4]Jl"9| and their dynamics are slightly complicated. There are some studies on the 
interaction between a cosmic string and a black hole with the hope of detecting gravitational waves from cosmic 
bJ[), strings. De Villier and Frolov studied [|]|| the dynamics of the scattering and capturing process of an infinitely thin 
test string by a Schwarzschild black hole. However less is known on the interaction between a domain wall and a 
black hole; only the static configurations have been studied so far. Morisawa et al. [^6| showed the existence of static 
configurations of a thick domain wall intersecting a Schwarzschild black hole in its equatorial plane. Christensen 
et al. || numerically found a series of static configurations of a domain wall as a Nambu-Goto membrane in the 
Schwarzschild spacetime, some of which represent an intersecting pair of a domain wall and a black hole. When 
arranged in sequence, these configurations seem to represent a scattering and capturing process of a domain wall by a 
black hole. In most situations of cosmological interest, the thickness of a domain wall will be much smaller than the 
horizon radius of a solar-mass or primordial black hole, so that the infinitely thin wall (membrane) approximation 
will be valid in such cases. When we further neglect the gravitational effect of a domain wall, its spacetime history, 
i.e., the world sheet is governed by the Nambu-Goto action and it is called the Nambu-Goto membrane. 

A membrane lying in the equatorial plane of the black hole spacetime has the highest symmetry among the config- 
urations which represent an intersecting pair of a membrane and a black hole. This simple configuration is a possible 
candidate for a final state of the scattering and capturing process provided such a state stably exists. Having this 
in mind, we study the stability of a Nambu-Goto membrane at the equatorial plane of the Reissner-Nordstr0m-de 
Sitter (RNdS) spacetime. The charge of the black hole and the positive cosmological constant are included to see 
their effect on the stability. This mimics the situation where a domain wall resulting from a vacuum phase transition 
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is intersecting a charged primordial black hole in the early universe. We consider the linear perturbation of the 
membrane by means of the covariant perturbation formalism J^Jl^|. The perturbation equation reduces to a wave 
equation on the world sheet with the mass term, which is negative in the presence of the charge of the black hole 
or the positive cosmological constant. Accordingly instability is naively expected. Solving the perturbation equation 
reduces to solving a two point boundary value problem. We numerically solve it by the shooting and the relaxation 
methods. We find that a membrane intersecting a RNdS black hole is unstable. We find the instability at least in 
the range of 0.001 < e/M < 1 for A > 0, and 0.75 < e/M < 1 for A — 0, whereas the membrane in the Schwarzschild 
background is stable. 

While our treatment of a domain wall is simple with both its thickness and gravitational effect neglected, our 
model of the membrane-black hole system is all the more elementary and relevant in other contexts. The interaction 
between black holes and extended objects is important as an elementary physical process. Indeed much attention has 
recently been paid to the interaction between black holes and extended objects such as membranes or strings in the 
string/M-theory ]9|-|l2]], the brane world scenario p],p|,[7|,p7 18 etc.. For instance in the context of the brane world 



models, our universe is described as a four-dimensional domain wall in the five dimensional bulk spacetime. Some 
people argue that a black hole on the gravitating membrane is realized as a "black cigar" in the bulk spacetime which 
intersects the membrane [Q . 

The rest of the paper is organized as follows. In Secjl| we review the covariant perturbation formalism and specify 



the perturbation equation and the boundary conditions. The numerical algorithm is also explained there. In Sec. Ill 
we present the results of numerical and analytical considerations. Finally we summarize and discuss the results in 
Sec" 



IV 



II. BASIC EQUATIONS 
A. Equation of motion 

The history of a membrane (world sheet) is described by the timelike hypersurface (S,7 a fc) embedded in the four- 
dimensional spacetime (M, c/^,,). The embedding is given by 

x»=XHj?) M = 0,---,3, a = 0,---,2, (1) 

where £ M 's are the spacetime coordinates and £ a 's the world sheet coordinates. The induced metric j a b on £ is given 
by 

lab = X%X» b9llv . (2) 
The dynamics of the membrane is described by the Nambu-Goto action 

S[X",X»] = -< r f d 3 ^, (3) 

where a represents the surface energy density of the membrane in its rest frame. 
We consider the variation of the action with respect to X^, 

' )S a(T^ b X%X^ b + axn, (4) 



sx. 



where T^o is the spacetime affine connection and □ is the d'Alembertian on the world sheet 



□ = 



1 d , ab d 



-ae^W ] - (5) 



Recall the Gauss- Weingarten equation 



D *>X% + K X % X fi = K ^ , (6) 



where K ab is the extrinsic curvature defined as 

K ab = -MV„n„ (7) 
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n M is the normal to the world sheet, D a is the world sheet covariant derivative and V M is the spacetime covariant 
derivative. Contracting Eq. (||) with 7 Qb , we obtain 

nX» + T^ ab X%XP = ^Katn^ . (8) 

With this equation Eq. (||) becomes 

" S a 7 ab K ab n^ = 0. (9) 



5X^ 



Eq. (||) has only the component perpendicular to the world sheet. Hence the variation parallel to the world sheet has 
no physical meaning. Finally the equation of motion of the membrane is 

K = "f ab K ab = 0. (10) 

This is the equation of minimal surfaces. In general, Eq. ([l0]) cannot be solved analytically. 

B. Perturbation equation 

As noted above, the physically meaningful measure of the perturbation SX^ is the scalar 

$ = n^SX*. (11) 
Taking the variation of Eq. ([To]) with respect to $ and setting it to zero, we obtain the perturbation equation 

□ $ - ( 3 R - A R 11V W 1V ) $ = 0, (12) 
where = g M „ — n M n„ is the projection operator onto the world sheet. 

C. Perturbation of a membrane at the equatorial plane of the RNdS spacetime 

A membrane lying in the equatorial plane of the black hole spacetime has the highest symmetry among the config- 
urations found by |3j which represent an intersecting pair of a membrane and a Schwarzschild black hole. This simple 
configuration is a possible candidate for a final state of the scattering and capturing process of a membrane and a 
black hole provided such a state stably exists. 

In general, spherically symmetric black hole spacetimes have the discrete symmetry about the equatorial plane, 
i.e., the equatorial plane consists of fixed points of this symmetry. Hence the equatorial plane is totally geodesic and 
K = 0. 

The most general spherically symmetric black hole spacetime is the Reissner-Nordstr0m-de Sitter spacetime 

g^dx^dx" = -fdi 2 + f^dr 2 + r 2 d6 2 + r 2 sui 2 6dip 2 , (13) 

where 

2M e 2 A 2 

/ = 1 + - " T-r 2 , 14 

r H 3 

A is the cosmological constant and M and e are the mass and charge of the black hole, respectively. The relation 
between the spacetime coordinates {t, r, 9, (p} and the world sheet coordinates {T, R, of the membrane at the 
equatorial plane is simply 

t = T, r = R, 0=|, ip = V. (15) 

Then the induced metric (|2|) is 

JabdCd£ b = -fdT 2 + f^dR 2 + R 2 d^> 2 . (16) 

This means that the geometry of the membrane is a three-dimensional black hole spacetime. The perturbation 
equation (n2h becomes 
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A $ = 0. 



(17) 



The mass term becomes negative when the charge of the black hole or the positive cosmological constant is present, 
which implies instability. This is why we consider the equatorial plane of the Reissner-Nordstr0m-de Sitter spacetime. 
In addition, it mimics a cosmologically interesting situation where a domain wall resulting from a vacuum phase 
transition is intersecting a charged primordial black hole in the early universe. 

By a separation of variables, $ = (x(r)/y / r) exp(zo;t + irmp) , Eq. ( |l7j ) is transformed into a Schrodinger type 
equation in the three-dimensional black hole spacetime 



d X 
dri 



V(r) X 



(18) 



V(r) = f 



m 2 e 2 \ f-2fr 



4r 2 



where r* = J r dr/f is the usual tortoise coordinate. For the cases M/0 we normalize r, e, A, u> as 

r/M -> r, e/M -> e, AM 2 -> A, ujM -> w. 



(19) 



(20) 



D. Boundary conditions 



The present problem reduces to solving a Schrodinger type equation, and so does the problem of the metric 
perturbation of a four-dimensional black hole spacetime. Hence we will follow the stability analysis of the latter. 

We consider the cases M ^ 0. Then the tortoise coordinate r* goes from — oo to +oo. r* = — oo corresponds to the 
event horizon and r* = +oo to infinity (A = 0) or the cosmological horizon (A > 0). The potential (fL9|) approaches 
zero as r* goes to ±oo. Therefore the asymptotic solution to Eq. ( fl8| ) is a linear combination of e +luJT ' and e~' lujr * . 
With the time dependence of the solution of the form e lut , the solutions e +lur ' and e~~ zlJjr * represent an ingoing wave 
and an outgoing wave, respectively. At infinity or the cosmological horizon we admit an outgoing wave, and at the 
event horizon an ingoing wave. This is because we assume no external source of perturbation of the membrane during 
its evolution. Hence we set the boundary conditions as follows 



X 



C' ou t exp(-iwr* 
C in exp(-Hwr») 



(r* 
(n 



-co) 
-oo) 



(21) 



with Ci n i out constant. This form of boundary conditions appears in the standard analysis of a quasi- normal mode of 
a black hole [g]. 

With the time dependence e Mt , Im(u>) < corresponds to unstable modes. Then from the boundary conditions 
(pl|), an unstable solution to Eq. ( fr§| ) decays as e,-\ Im (. u )\ r * when r* goes to +oo and as e' /m ^" r * when r* goes to 
— oo. Multiplying Eq. ( [18] ) by \ (upper bar means the complex conjugate) and integrating by parts, we obtain 



dx 



dr* 



+ V\ X \ dr* + 



, dx 



UJ 



-t-oo 



Ixl dr*. 



(22) 



As long as we consuder unstable modes, the surface term on the l.h.s. 
converge. Therefore we obtain 



vanishes and the integrals on both sides 



LU 2 = 



dx 
dr* 




\x\ 2 dr* 



(23) 



Since the r.h.s. of Eq. (|23j) is real, to is real or pure imaginary. Here we seek for unstable modes (Im(u>) < 0). So we 
can set to — ia (a < 0). Then the boundary conditions ( pi|) read 



Cout exp(+crr*), (r* — > +oo) 
Ci n exp(-crr*), (r* — ► -oo) 



(24) 



We shall examine the eigenvalue problem of Eq. (08) subject to the boundary conditions (p4|). 
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E. Algorithm 



The eigenvalue problem considered is reduced to a two point boundary value problem. We define four dependent 
variables as 

^X 

2/i = X, 2/2 = -j— , 2/3 = cr, 2/4 = A, (25) 

where all j/i's are functions of r*, and cr and A are constants. A is the ratio of C;„ to C out (or C OMt to C% n ). The 
evolution equations are 

y[ = 2/2, 2/2 = [V + (1/3) 2/3 = 0) v'i = 0, ( 26 ) 

where the prime denotes the differentiation with respect to r*. We impose two sets of boundary conditions. The set 
1 (2/4 = C out /C m = Ax) is 

2/1 = cxp(~y 3 r* 1 ) ,y 2 = -j/ 3 2/i , (at r* = r*i) (27) 
yi = y 4 exp(y 3 r* 2 ) , y 2 = 2/32/1 (at r* = r* 2 ). (28) 

The set 2 (y 4 = C in /C out = A 2 ) is 

2/1 = 2/4exp(-y 3 r» 1 ) , y 2 = -2/32/1 i (at r* = r* x ) (29) 

2/1 = exp(y 3 r» 2 ) , y 2 = y 3 yi (at r* = r* 2 ). (30) 

The condition A\ x A 2 = 1 is used for checking the reliability of calculations. Note that the boundary conditions 
( p7[ - |30| ) leave ?/ 3 , j/4 arbitrary. In general an arbitrary choice of y 3 , j/4 at one boundary and the subsequent integration 
of the dependent variables do not ensure that the boundary conditions at the other boundary are satisfied. We solve 
this problem by the shooting method. We also perform the calculation by the relaxation method to confirm the results 
of the shooting method. 



III. RESULTS 
A. Schwarzschild case 



For pure Reissner-Nordstr0m cases the potential (|T^) is positive semi-definite for m > 1: 



V{r) 



1 



Air 2 



2r 3 4r 4 



> 



(r > r+), 



(31) 



where r + = l + ^/l — e 2 is the radius of the outer event horizon. The equality is satisfied only if m = e = 1. Therefore, 
by the argument of the energy integral 0, there is no unstable mode subject to the boundary conditions (|24| ) when 
m > 1. 

For the Schwarzschild case with m = the radial function R(r) = x( r )/V^ satisfies 



1 - 



2\ d 2 R ldR 



dr 2 r dr 1 — 2/r 

In terms of a new variable £ = ln(r — 2) this equation becomes 

d 2 R 9 9 _ 
- — + * 2 r 2 R = 0. 



R = 0. 



(32) 



(33) 



The solution to this equation is regarded as the zero energy eigenfunction for the potential a 2 r 2 . This potential is 
positive definite since cr is real and negative. Therefore we have no solution which decays when £ approaches ±00 
(the event horizon and infinity) as required by the boundary conditions (^4|) 



R(r) = 



X(r) 



(r — > 00 — > +00)) 
(r-2)- 2CT /V2, (r ^ r+ (C ^ -00)) 



(34) 



Hence there exists no unstable mode for the Schwarzschild case. 
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B. RNdS case 



Before describing the numerical results, we comment on the upper bound of to for which unstable modes may exisit 
in the Reissner-Nordstr0m-de Sitter case. The potential (|l9|) takes the form 



where / is defined by Eq. (jbj 



and 



v(r ,{fF{r)l{^) (e = 0) 

V( >~ \/G(r)/(4r 4 ) (e ? 0) ' ^ 



F(r) = -5Ar 3 + (4m 2 - l)r + 6 (36) 



G(r) = -5Ar 4 + (4m 2 - l)r 2 + 6r - 9e 2 . (37) 



Firstly we consider e = case. F(r) is a monotonously decreasing function of r and has the only one root for r > 0. 
Therefore, by the arguement of the energy integral Q , the sufficient condition for the stability of perturbation is 

F{r c ) > 0, (38) 

where r c is the radius of the cosmological horizon. Eq. (pq) reads 



2 r c - 6 + 5Ar e 3 9 

to > =4 , (39) 

4r c r c 

where /(r c ) = is used. When e = 0, r c > 3. Therefore Eq. ( |39| ) shows that at least for to > 2 the potential ( |l9| ) is 
positive definite and that a membrane is stable against such perturbation. 

Next we consider e ^ case. For r > 0, G(r) either stays negative, or becomes positive for an interval {r\,r2) but 
otherwise negative. Then the sufficient condition for the stability of perturbation is 

G(r+) > and G(r c ) > 0, (40) 

which excludes the possibility that G{r) stays negative, and ensures r± < r+ < r c < ri- When e 2 increases, r+ and 
r2 decrease and T\ and r c increase. Hence it is sufficient to consider e 2 — 1 case. In this case Eq. (^) reads 

, f 16r c - 3r c + 24 16r+ - 3r+ + 24 1 f 3 6 3 6 1 

m 2 >max c — -£ , + - 2 + } = 4 + max^ -— + -— + }, 41 

[ 4r c 2 4r + 2 J [ 4r c r c 2 4r+ r+ 2 J 



where again f(r+) = f(r c ) = are used. When A > 0, r + > 1 + yl — e 2 = 1 (now e 2 = 1). For r > 1 

-t|j<-!- + -!< 5 - 25 - (42) 
128 4r r 2 

Therefore Eq. ( |4l| ) shows that at least for to > 4 there exists no unstable mode. 

The condition that the Reissner-Nordstr0m-de Sitter spacetime has a static region is A < A max , where 



A_ (e 2 ) = 3 ( 1 + ^ZE1 . (43) 

V ; (18-4e 2 )(3 + V9-8e 2 )-24e 2 V ; 

This is a monotonously increasing function of e 2 and has the minimum value A max (0) = 1/9. So we numerically 
investigate the cases AM 2 = 0, 1/1000, 1/100, 1/10. We found unstable modes for various parameter sets (A, |e|). The 
eigenvalues of to = modes are shown in Fig. Vu. The results of the shooting method and the relaxation method agree 
very well, the error being within 5%. 

Fig. [l] shows that when |e| becomes larger, we have a larger \a\ and hence a higher growth rate of the perturbations. 
Therefore the charge of the black hole indeed destabilizes the membrane at the equatorial plane. Since the membrane 
is electrically neutral, this instability should be understood as the curvature effect of the charge of the black hole. 

As to the effect of the cosmological constant A, we find no simple tendency: for a given |e|, a larger A does not 
necessarily lead to a larger |er|. However with a non- vanishing cosmological constant the magnitude of the eigenvalue 
is typically 10~ 2 for 0.001 < \e\/M < 1 in contrast to the RN cases, in which the magnitude of the eigenvalue almost 



G 



exponentially decreases with decreasing specific charge and is far below 10~ 2 for |e|/M < 0.80. Therefore we can 
at least say that the cosmological constant also destabilizes the membranes. This expectation is supported by the 
consideration of the de Sitter background case (see Appendix [A|) . 

We plotted some of the eigenf unctions for A = cases in Fig. ||. Some features of the profiles are listed in Table |. 
When we decrease the value of the specific charge of the black hole, the peak of the corresponding profile gets higher 
and approaches the event horizon. 

As noted above we found that the charge of the black hole destabilizes the membrane in general. However, in the 
present analysis, we could not show the existence of unstable modes for small |e| (< 0.75) when A = due to the 
difficulty in the numerical calculation as follows. As |e| becomes smaller, \a\ decreases approximately in powers of 
|e| for 0.75 < |e| < 1 (See Fig. ||). Judging by the extrapolation of this relation to the regime |e| < 0.75, \a\ will be 
infinitesimal when |e| ~ 0. In order for the asymptotic forms e ±lujr ' of \ to be justified, |V(7")| must be much smaller 
than |o; 2 | = \a 2 \ at the boundaries. Then the coefficients in the evolution equations ( |26| ) get extremely small and we 
suffer from underflows of numerical calculations. By the same numerical difficulty, we leave open the possibility that 
unstable modes with m > 1 exist in the presence of the positive cosmological constant. 



IV. DISCUSSION 



In this paper we numerically studied the stability of a Nambu-Goto membrane at the equatorial plane of the 
Reissner-Nordstr0m-de Sitter spacetime and found that in general such a membrane is unstable when the black hole 
is charged. 

A membrane is a two-dimensional extended object in the three-space. So when it moves towards a black hole, it 
is expected to be inevitably captured by the black hole. A series of configurations found by Christensen et al. seem 
to represent such situations. Among these configurations a membrane lying in the equatorial plane of the black hole 
spacetime has the highest symmetry and is most likely to be the final state of the scattering and capturing process. 
On seeing our result, however, we expect that it is not the case. In particular, Fig. |^ allows us to imagine that 
the membrane moves away from the equatorial plane. Then what eventually happens to that membrane? Causality 
prohibits the membrane from escaping from the event horizon. By analogy with the scattering process of a cosmic 
string off a black hole [^,^) the membrane might experience large deformation and the topology of the membrane might 
change. Here are some speculations on the fate of the membrane; it may settle down to some other configuration than 
the equatorial plane, or it may break up into two parts with one swallowed by the black hole and the other escaping 
to infinity. We need to perform a full dynamical computation to resolve this issue. However that is beyond the scope 
of this paper. 

We ignored the gravitational effect of a domain wall whereas a gravitating domain wall is known to make a repulsive 
gravitational field |Q, which is opposite to the strong attractive gravity of a black hole. It is quite intriguing to find 
out the consequences of the competition of these opposite forces. 
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APPENDIX A: DE SITTER AND MINKOWSKI BACKGROUND CASES 



To understand the destabilizing effect of the positive cosmological constant, we compare the de Sitter and Minkowski 
cases. In these cases Eq. ( |l8| ) is rewritten in terms of a new function P(r*) = xi r *) I \fr*- a s 

d 2 P 1 dP ( 1 1 7 A 2 5 ., 2 \ _ „ 

= h — 7T ~ A + fTH AV ) P = 0. (Al) 

dr 2 r*dr* \4r 2 4r 2 6 12 J y ' 

where m is set to zero for simplicity. The reason we use P{r*) is that in the de Sitter and Minkowski cases the 
potential (|l9|) for x( r *) diverges at the center r = and that numerical calculations become unstable. 
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In the Minkowski case = r and Eq. (]AlJ) becomes Bessel's differential equation of order 

d 2 P 1 dP 



d(ar) err d(ar) 



P = 0, (A2) 



the solution of which is a linear combination of the modified Bessel functions I^iar) and Ko(ar). However Jo (or) and 
Ko(crr) diverge at r — > oo and r = 0, respectively and do not satisfy the regularity conditions. Hence membranes are 
stable in the Minkowski case. (The stability for m > 1 cases is also verified.) 

When the positive cosmological constant is present, the situation changes. In the limit of r — > (7-* — > —00), 



Eq. ( [Al| ) reduces to 

d 2 P 1 dP /ll 



(//•- ' r* '//•. V 9 



Unlike Eq. (A2), the solution to Eq. (fO) differs depending on the sign of (11 A/9 - a 2 ). When \a\ < ^/llA/9, 



the solution is a linear combination of the Bessel function Jo (or*) and the Neumann function Ao(ar*), where a 
\/|HA/9 — cr 2 |. When |cr| > v/llA/9, the solution is a linear combination of the modified Bessel functions Jo(ar*) 
and Jfo(aT*). The regularity condition at the center excludes Nq and Kq as the asymptotic solutions. 

When \a\ < -y/llA/9, the asymptotic behavior of P near the center changes from the Minkowski case and we expect 
the emergence of unstable modes since Jo is a bounded function. In fact we found unstable modes by numerical 
calculations. The method is similar to the one for the RNdS cases. The asymptotic form of P near the cosmological 
horizon r* — * 00 is determined to be e CTr */y / rv by the boundary condition similar to Eq. (|^). The eigenvalues are 
-0.18257, -0.05773, -0.01826 for A = 1/10, 1/100, 1/1000, respectively. Though we performed numerical calculations 
just for three values of the cosmological constant, unstable modes are expected as long as the positive cosmological 
constant is present. 



APPENDIX B: CODE CHECK 



We calculated the energy spectra for a one-dimensional square- well potential to check the reliability of our numerical 
code. 

We consider an eigenvalue problem 

-^<l> + V{x)<t> = -o 2 <l> (Bl) 

for a square-well potential 

v(x) = \- 1 (N<^>o) . (B2) 

v ' 10 (otherwise) v ' 

The problem reduces to a two point boundary value problem at the boundaries x = x\,X2 {x\ < —w,%2 > w). The 
dependent variables are 

d(f> 

Vi = <P, V2 = -j-, 2/3 = 0-, t/ 4 = A, (B3) 
dx 

where a and A are constants. The evolution equations and the boundary conditions are similar to Eq. (p6| ) and Eqs. 
©H) or Eqs. (HIH), respectively. 
We can analytically show that the eigenvalues for the potential ([B^) satisfy the equation 

wr-2i<;s-2arcsms = (n = 0, ±1, ±2, • • •) (B4) 



where s = yl — a 2 . The eigenvalues a wn of Eq. (Bl) are labeled by w and n. We found good agreements of eigenva lue s 
a wn obtained by solving the two point boundary value problem, with those obtained by directly solving Eq. (B4) 
(Table |J. 
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APPENDIX C: CODE CHECK 2 



In this paper we present the results for which the shooting and relaxation methods agree within 5 % error. The 
agreement is, however, far better than 5% in most cases. The comparison of the two methods is summarized in 
Table III for some cases. 
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FIG. 1. The plots of the eigenfrequency, |crM|, versus the specific charge, \e\/M, of the black hole for the m = cases. The 
top left panel shows the plot for the AM 2 = case, the top right for the AM 2 = 1/1000, the bottom left for the AM 2 = 1/100 
and the bottom right for the AM 2 — 1/10. We obtained two sequences of the eigenvalues except for the AM 2 = case. 
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FIG. 2. The radial profiles R(r) = x{ r )/Vr °f the m = mode eigenfunctions for the pure Reissner-Nordstr0m cases 
(\e\/M — 1, 0.99, 0.90, 0.80, 0.75). The abscissa is normalized in the unit of the horizon radii r+. 
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TABLE I. Some features of the radial profiles in Fig. 2. 



charge (|e|/M) 


event horizon (r+/M) 


peak location r pea ^/M 


peak amplitude -R(r pae k) 


r peak/ r + 


1 


1 


1.202 


0.22 


1.202 


0.99 


1.141 


1.281 


0.57 


1.123 


0.90 


1.436 


1.468 


0.79 


1.022 


0.80 


8/5 


1.603 


0.79 


1.002 


0.75 


1.661 


1.662 


0.78 


1.001 



0.1 



0.01 



6 0.001 

c 



0.0001 




1e-05 1 1 1 1 

0.75 0.8 0.9 1 

In |e/M| 

FIG. 3. The plots of the eigenfrequency, ln|crM|, versus the specific charge , ln|e/M|, of the black hole for the AM 2 = 0, m — 
cases. We find an approximate linear relation between In \crM\ and In \ e/M\. 



TABLE II. Code check. Comparison of numerically obtained eigenvalues for the square-well potential (B2) with those 
obtained from Eq. (B4). 





analytic (from Eq. (pll|) ) 


numerical 


A (see Eq. (p^)) 


parity 


on 


0.673612 


0.67361 


+1 


even 


(T21 


0.85723 


0.85723 


+1 


even 


C22 


0.319023 


0.31902 


-1 


odd 


0-31 


0.920798 


0.92080 


+1 


even 


C32 


0.650366 


0.65037 


-1 


odd 


(741 


0.949724 


0.94972 


+1 


even 


C"42 


0.785671 


0.78566 


-1 


odd 


C43 


0.438306 


0.43830 


+1 


even 


C51 


0.965261 


0.96527 


+1 


even 


CT52 


0.854684 


0.85467 


-1 


odd 


C53 


0.641057 


0.64105 


+1 


even 


CT54 


0.192693 


0.19269 


-1 


odd 
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TABLE III. Comparison of the shooting and relaxation methods. Some selected eigenvalues oM are listed. 



Specific charge (|e|/M) 


Shooting 


Relaxation 


AM" = 




1 


O.oOl X 1U 


f\ ^98 v i n~ 2 

U.OZO X 1U 


0.99 


A (KAA v 1 fl — 2 


A f\AC\ v i n~ 2 
4.04U X 1U 


0.97 


9 Q98 v 1 n — 2 
z.yzo x iu 


9 Q97 v 1 0~ 2 


0.95 


1 QOO w i n— 2 


i Q9Q w i n — 2 
l.yzo X 1U 


0.93 


1 9£8 v i n — 2 


i 9fiK v i n — 2 


0.91 


8 1 Qfi x in -3 


8 182 x 10~ 3 


0.89 


5.139 x 10~ 3 


5.120 x 10~ 3 


0.87 


3.104 x 10~ 3 


3.084 x 10~ 3 


0.85 


1.797 x 10~ 3 


1.778 x 10~ 3 


0.83 


9.939 x 10~ 4 


9.792 x 10~ 4 


0.81 


5.242 x 10~ 4 


5.180 x 10~ 4 


0.79 


2.632 x 10~ 4 


2.559 x 10~ 4 


0.77 


1.255 x 10~ 4 


1.211 x 10~ 4 


0.75 


5.676 x 10 -5 


5.416 x 10~ 5 


AM 2 = 0.1 




1 


9.9727 x 


10"^ 


9.9764 x 10^ 


l(2nd seq.) 


7.2435 x 


10~ 3 


7.3476 x 10~ 3 


0.99 


9.6971 x 


10~ 2 


9.7009 x 10~ 2 


0.90 


8.1253 x 


10~ 2 


8.1291 x 10~ 2 


0.80 


7.0378 x 


10~ 2 


7.0409 x 10~ 2 


0.70 


6.2124 x 


io- 2 


6.2149 x 10~ 2 


0.60 


5.5387 x 


10~ 2 


5.5406 x 10 -2 


0.50 


4.9761 x 


10~ 2 


4.9775 x 10~ 2 


0.40 


4.5092 x 


10~ 2 


4.5103 x 10~ 2 


0.30 


4.1352 x 


10~ 2 


4.1361 x 10~ 2 


0.20 


3.8587 x 


10~ 2 


3.8594 x 10~ 2 


0.10 


3.6878 x 


10~ 2 


3.6885 x 10~ 2 


0.01 


3.6305 x 


10~ 2 


3.6311 x 10~ 2 


0.001 


3.6298 x 


10~ 2 


3.6305 x 10~ 2 
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